511 KeV Photons From Color Superconducting Dark Matter 
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We discuss the possibility that the recent detection of 511 keV 7 rays from the galactic bulge, 
as observed by INTEGRAL, can be naturally explained by the supermassive very dense droplets 
(strangelets) of dark matter. These droplets are assumed to be made of ordinary light quarks 
(or antiquarks) condensed in non-hadronic color superconducting phase. The droplets can carry 
electrons (or positrons) in the bulk or/and on the surface. The e^e~ annihilation events take place 
due to the collisions of electrons from the visible matter with positrons from dark matter droplets 
which may result in the bright 511 KeV 7-ray line from the bulge of the Galaxy. 

PACS numbers: 98.80.Cq, 95.30.Cq, 95.35.+d, 12.38.-t 



Introduction. — The recent detection by the SPI spec- 
trometer on the International Gamma-Ray Astrophysics 
Laboratory (INTEGRAL) satellite of a bright 511 KeV 
7-ray line from the bulge of the Galaxy with spherically 
symmetric distribution []| has stirred the research of the 
fundamental physics that describes the cosmological dark 
matter. 

The flux of 511 KeV photons (with a width of about 
3 KeV), produced by thermalized electron positron pair 
annihilation processes, has been measured to be 9.9^2 1 ^ 
10^** photons cm^^ s~^ and has an angular distribution 
with half maximum at 9° (6° to 18° at 2cr confidence), in 
good agreement with previous measurements 2] . 

The source of these thermalized positrons in the bulge 
of the Galaxy have been the subject of much debate. 
Some proposals suggest astrophysical processes, includ- 
ing neutron stars or black holes [3 , pulsars |j| , radiactive 
nuclei from supernova |5J or cosmic ray interactions with 
the interstellar medium |g], but it is rather uncertain 
which fraction of positrons produced in such processes 
can escape and, moreover, how they could fill the whole 
galactic bulge J^. 

Recently it has been discussed that light dark mat- 
ter particles annihilating into e"'"e^ pairs in the galactic 
bulge may be the source of the thermalized positrons that 
produce the 511 KeV emission line Q, see also related 
works, [2|, [ill, [ll|. The shallow density distribution 
of dark matter in the bulge of the galaxy p{r) '^ ?" '', 
with 7 = 0.4 to 0.8, explains very naturally the angular 
distribution of detected 511 KeV 7 photons. 

Dark matter as color superconductor. — We want to 
elaborate the proposal [g] in the context of a cosmolog- 
ical scenario when dark matter consists of very dense 
(few times the nuclear density) macroscopic droplets 
of ordinary light quarks ( or/ and antiquarks [13, [l^ ) 
condensed in non-hadronic color superconducting phase, 
similar to the Witten's strangelets [lj|. 

In this Letter we will argue that color superconducting 
dark matter also provides a natural and simple frame- 
work to explain the detected emission of 511 KeV pho- 
tons from the galactic bulge with the appropriate angular 
distribution and intensity. Indeed, the main required in- 



gredients of the proposal are automatically present in our 
scenario: a large number of positrons is always present in 
antimatter dark matter droplets, see below. We argued 
in 13] that chunks of quarks or antiquarks in condensed 
color superconducting phase may be formed during the 
QCD phase transition and they may serve as dark mat- 
ter. This scenario is based on the idea that while the 
universe is globally symmetric, the antibaryon charge can 
be stored in chunks of dense color superconducting (CS) 
antimatter. In different words, the baryon asymmetry of 
the universe may not necessarily be expressed as a net 
baryon number if the anti-baryon charge is accumulated 
in form of the diquark condensate in CS phase, rather 
than in form of free anti-baryons in hadronic phase. We 
explained in 13' why such a scenario does not contradict 
the current observational data on antimatter in the Uni- 
verse. This is mainly due to the very small volume occu- 
pied by dense droplets and specific features of interaction 
between color superconducting phase and conventional 
hadronic matter. We also argued that the observed cos- 
mological ratio between the energy densities of dark and 
baryonic matter, floM '^ ^b within an order of magni- 
tude, finds its natural explanation in this scenario: both 
contributions to fl originated from the same physics at 
the same instant during the QCD phase transition. As 
is known, the relation Qb ~ ^dm between the two very 
different contributions to Q is extremely difficult to ex- 
plain in models that invoke a DM candidates not related 
to the ordinary quark/baryon degrees of freedom. The 
baryon to entropy ratio UB/n^ ~ 10~^° would also be a 
natural outcome in this scenario. We refer to the original 
papers [13I, [13 for the details. Here we want to mention 
only the fact that the baryon charge of massive droplets 
does not change the nucleosynthesis calculations because 
in the color superconducting phase it is not available for 
nuclearsynthesis when the baryon charge is locked in the 
coherent superposition of Cooper pairs. Therefore, while 
the massive droplets carry a large baryon charge, they 
do not contribute to Sis, but rather, they do contribute 
to the "non-baryonic" cold dark matter ^dm of the uni- 
verse [13, [ii- 

Before we estimate the probability of the e'^e~ anni- 



hilation which results in 511 KeV Hne, we would like to 
make a short review on basic properties of dense droplets 
in color superconducting phase, which will be referred as 
QCD balls [l^,[l3j in what follows. The color supercon- 
ducting state of quark matter is a novel phase in QCD 
that is realized when light quarks are squeezed to a den- 
sity which is a few times the nuclear density. The ground 
state in this phase is a single coherent state with diquark 
condensation, analogous to Cooper pairs of electrons in 
BCS theory of ordinary superconductors. In the approx- 
imation of three light quarks 77i„, mrf,ms <C n and rela- 
tively large chemical potential jjl ^ Aqcd, the so-called 
color-flavour-locking (CFL) phase is a preferred state of 
matter, see original papers jl^ and recent review IQj on 
the subject. For physical value of rus and fi ~ 500 MeV 
a number of different CS phases may result. It is not 
the goal of this letter to describe a variety of possibilities 
when parameters (such as nis and /i) vary. Rather, we 
want to emphasize below that the sufficient number of 
positrons will always accompany the QCD balls made of 
antimatter (QCD anti-balls). 

Indeed, first of all, consider the most symmetric, the 
CFL phase. While this phase does not support the lep- 
tons in the bulk 17,], the finite volume effects lead to the 
accumulation of the positive charge on the surface [ig, 
which must be neutralized by negative electron charge 
(for droplets made of matter) . For droplets made of an- 
timatter, the corresponding positron charge will be accu- 
mulated. In most other phases which may be realized in 
nature, the Icptons will be present on the surface as well 
as in the bulk of a droplet. The electron density can be 

roughly estimated as rie — j^ , with /Ue being the electron 
chemical potential (in case of matter droplets) or positron 
chemical potential in case of anti-matter droplets. In this 
case, the electrons (positrons) in droplets can be treated 
as fermi liquid. A numerical estimation of ^e strongly 
depends on the specific details of CS phase under con- 
sideration, and varies from few MeV to tens MeV, [l9j- 

However, the important property which plays essen- 
tial role for the present work (and which is shared by 
all CS phases), is as follows. Consider an electron which 
hits the DM droplet (made of anti matter). What is 
the fate of this non relativistic electron? It can form a 
bound state (positronium with arbitrary quantum num- 
bers |n, I, m)) which eventually decays to two ~ bllKeV 
photons. It may also annihilate with energetic positron 
into two photons in non resonance manner with emit- 
ting 2 7's with a typical energy determined by /ie ( few 
MeV scale). However, the probability for the later an- 
nihilation is suppressed by small coupling constant a^, 
in comparison with the former process, when the proba- 
bility for the formation of positronium from two nonrel- 
ativistic particles e"*" and e~ could be order of one. In- 
deed, the probability for the positronium formation (as 



well as for its decay to free e~^e~ pair) if the system gets 
an instantaneous jolt (with relative momentum q = mv) 
is determined by the overlap of two wave functions ^ 
KV'outlV'm)!^ ~ I /e~''/"e*'''''d^rp where e"''/" represents 
a typical positronium wave function in state \n, I, m) with 
a ~ 10^*^ cm. Of course, this expression assumes the va- 
lidity of instantaneous perturbation theory when parame- 
ter qa » 1, while the maximum probability is achieved 
when ga ~ 1, see below. It is obvious that the main 
contribution to the positronium formation is due to the 
process when the incoming electron picks up a positron 
from the droplet with a typical velocity determined by 
the condition: go ~ 1. This corresponds to v/c ~ a for 
a typical positronium size, a ^ h jme^ . Eventually, it 
decays to two ^ 511 KeV photons. The flux of emitted 
photons produced by this mechanism will naturally have 
a width of order T/{5llKeV) ~ v/c ~ a ~ 10"^, which 
is what observations apparently suggest [ij . We note that 
the positronium formation (with consequent emission of 
511 KeV photons) is expected to occur on the surface of 
the droplet such that considerable portion of 511 KeV 
photons leave the system without reabsorption. 

To conclude: the annihilation cross section for the elec- 
tron falling to the DM anti-droplet is given by the geo- 
metrical size of the object, A-kR^, while a typical width of 
outgoing flux of 511KeV photons is of order F ^ am ~ 
few KeV. These features are very universal, do not de- 
pend on specific details of the phase under consideration, 
and remain unaltered for all possible CS phases. With 
these remarks in mind, we estimate the e^e^ annihila- 
tion rate and the fiux of 511 KeV photons and compare 
it with the observational available data. 

First rough estimate — We start with a first estimation 
of the annihilation assuming that visible matter density 
follows the spatial distribution of dark matter, with the 
fixed ratio given by the cosmological ^b/^dm- We also 
assume that the electron density from the visible matter 
is roughly determined by the number density of protons. 
The system could be in ionized state (HII) or in neutral 
atomic hydrogen state. It is quite obvious that corre- 
sponding calculations would lead to strong underestima- 
tion of the annihilation rate because the visible matter is 
strongly peaked in the center of galaxy, the effect which 
is completely ignored in our first estimate. The positive 
elements of such an assumptions are: a) it allows us to 
follow closely the original analysis in 8], such that the 
spatial integration over matter density can be extracted 
from this paper, and the corresponding comparison with 
[8| can be made; b) it gives us a lower bound for the 
corresponding annihilation rate as argued above. More 
importantly, this lower bound depends only on a typical 
size of the droplets, and does not depend on specific as- 
sumptions on behavior of visible matter density in the 
center of galaxy. 

The estimation of the fiux of 511 KeV photons coming 
to Earth from the bulge of the Galaxy along the angular 



direction fi goes as follows. As we mentioned above, the 
number of electrons is roughly determined by the number 
density of protons, rif,- ~ n^, and all electrons which hit 
the QCD anti-ball (antidroplet made of antimatter) with 
radius R will annihilate such that a considerable portion 
of the process will lead to the production of two 511 KeV 
photons. The probability per unit time -J|- that this 
happens in the presence of a single QCD ball is given by 
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where vjc ^ 10~^ and we express the baryon density 
in terms of dark matter density, IGeV ■ ng c::^ pg c^ 
^b/^dmPdm — 0.15/9£ij\f to make our first rough esti- 
mate. In order to estimate the probability of such events 
per unit volume per unit time Jy, one should multiply 
eq.(0 by the inverse volume occupied by a typical QCD 
ball with a typical baryon charge B. In our framework 
when the dark matter is identified with QCD balls and 
anti balls with typical mass M ~ niN ■ B, the corre- 
sponding number density of the DM particles is nothing 
but riDM — iGeV 'E ■ Therefore, we arrive to the following 
estimate, 



our assumption on visible matter density distribution, 
Pb — O.lBpDM with pdm ~ r~'^ and 7 = 0.6, normalized 
to the local density pdm — Q-iGeV/cm? would lead to 
the total visible material (within 8.5 kpc region ) of about 
4 • IO^Mq instead of observed ~ IO^Mq. 

Nevertheless this simple estimate is very instructive. 
First of all, one can explicitly compare our expression 
(PJ with the corresponding formula from ref. *p\ when 
the same factor describing DM distribution, J(Ail)Ai7 
enters the relevant formulae. Secondly, even the obvi- 
ously underestimated expression (0 is not in contra- 
diction with the existing bound on such kind of dense 
droplets, see eq. (20) in ref. l?] where bound B > 10^° 
is quoted. 

Unaccounted effects. Further complications — Here we 
want to discuss some new effects (ignored above) which 
certainly increase the rate. Unfortunately, the corre- 
sponding estimates are strongly model dependent, see 
below, and, therefore, should be taken with some cau- 
tious. First, let us take into account the properties of 
the visible matter distribution in the galaxy in a more 
appropriate way than it is done above. We replace for- 
mula ^ by the following expression. 
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The total flux of photons resulting from annihilation is 
obtained by integrating eq. l(2Jl over the line of sight and 
over the whole solid angle of observation. The numerical 
evaluation was done in 8]. We follow their analysis and 
implement it in our framework. We arrive at 

$ = I ds I dn^^ ~ 

Ln dVdt 
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where J(Arj)Afi ee /^^ dnj{n) with 



(3) 



The number density of electrons and the number density 
of dark matter particles are estimated as before, n^- ^ 
UB ^ (lify), noM ^ lEwh- We parameterize DM 
density as 
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normalized to the local density pdm — O-^iGeV/cm^, 
which is the central value adopted by |3] . For the visible 
matter we adopt the following scaling behavior ( close to 
the r~^ behavior of an izothermal sphereJ22|), 



J(0) = ( 



0.3GeV/cm^' S.bkpc 



ds[pDM{s)]. (4) 
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In expression ((SJ we traded R from eq. lO in favor of 
B ~ 3 "n^cs assuming that a typical baryon number 
density in color superconducting phase, ncs, is three 
times the nuclear saturation density, ncs — 3no with 
no - (108Mey)3. 

The factor J(J1) has been evaluated in reference '0^ for 
different density profiles ^(r) ex r~^ with 7 = 0.4 — 0.8 
providing the best fit. For these favorite 7's the value 
J(AJ7)Ar2 has been shown to vary in the range 0.3 — 1.6. 
This value should be substituted into eq. © and 
compared with the observations, 9.9^^2! ^ 10^^ pho- 



normalized to the total visible mass of Mtot = 
J ' '^^ d^xpB — IO^^Mq within S.bkpc. We notice that 
such a peaked distribution of visible matter would, in 
principle, produce a narrower distribution of 511 KeV 
photons than currently preferred by observational val- 
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However, if we take 7 ~ for the dark matter, the an- 
gular distribution which follows from eq. H5I7|I would be 
close to the upper bound of the preferred value j3|. 

Combining eqs. H5I6I7(I we arrive to the following final 
expression for the flux 



As we mentioned above, we consider this estimate as 
the lowest extreme case (within our framework). Indeed, 
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In obtaining the estimate lO we cut off the integral 
/o 5 c '^^ *^ 0.5pc at small distances where the visible 
matter rises very fast ^ r~^-^ while DM behavior at such 
scales is absolutely unknown. Such a cut off obviously 
brings a large uncertainty into our estimate. There is 
also large uncertainty due to the unknown scaling prop- 
erties of the dark matter distribution at small distances. 
Finally, different clumps and structures (such as stars, 
MACHOS, astreroids, etc) of the baryonic matter can 
strongly enhance the estimate ^ due to the fact that 
a large number of positrons from the bulk (rather than 
from the surface) of the QCD balls can participate in 
annihilation. Unfortunately, we do not know how to 
account this effect properly. The main goal here is to 
demonstrate the sensitivity of the calculations with vari- 
ation of the visible and DM matter distributions: the 
difference between two estimates, Q and (jHl is almost 5 
orders of magnitude. 

Conclusion — The main goal of the present letter is to 
argue that the color superconducting dark matter ( intro- 
duced with quite different motivation !l2!|, 'iS') provides 
a natural and simple framework to explain the detected 
emission of 511 KeV photons from the galactic bulge with 
the appropriate angular distribution and intensity. While 
there are many other possibilities to explain this rate 
based on some specific DM features, such as annihila- 
tion or decay, the present proposal is unique in many 
respects and can be easily discriminated from other ex- 
planations based on DM particles. Indeed, a unique fea- 
ture of our scenario is proportionality of the local flux 
of photons to both the density of visible as well as dark 
matter, see eq.ljSJ. In other DM based explanations the 
local flux does depend only on the distribution of dark 
matter. The corresponding matter distributions are ob- 
viously very different. In particular, an observation of 
the effect on the same level but in a different direction 
(not pointing to the center of the galaxy) would rule out 
our explanation. 

We also point out that qq annihilation might be suffi- 
ciently large for relati vely energetic protons (with kinetic 
energy about IGeV) [lj|- In this case e+e^ annihilation 
with single bright 511KeV line (discussed in the present 
paper) would be accompanied by the wide (70 MeV -1 
GeV) 7 spectral density due to the baryon- antibaryon 
annihilation. These very different spectra in different fre- 
quency regions must be related to each other due to their 
common origin. Corresponding calculations are beyond 
the scope of the present work; however, a very simplified 
estimate of the corresponding flux can be obtained by re- 
placing electron velocity v in formula (jSJ by a proton ve- 
locity Vp/v ^ y/me/nip ~ 2 • 10^^ 23]. This corresponds 
to the assumption of the thermal equilibrium between 
electrons and protons in the ionized region in the bulge 
of the galaxy (the HII has a vertical scale hight of ^ 90 
pc |22|'). Estimated in such a way flux is definitely not in 



immediate contradiction with observations, where some 
access of 7 rays indeed has been observed by EGRET. 
We should add that the observed access has been inter- 
preted in 1241 ■ | 25l as due to the dark matter annihila- 
tion, and in l2_C| as due to pp annihilation. One more 
phenomenological consequence of the suggested scenario 
is that baryon- antibaryon annihilation which always ac- 
company bllKeV line eventually may be responsible for 
a "non observation" of the cusp behavior near the Galac- 
tic Center. It might be worthwhile to investigate this 
possibility in more details in future. 
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